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ABSTRACT 

The rotation rate in pre-supernova cores is an important ingredient which 
can profoundly affect the post-collapse evolution and associated energy release 
in supernovae and long gamma ray bursts (LGRBs). Previous work has focused 
on whether the specific angular momentum is above or below the critical value 
required for the creation of a centrifugally supported disk around a black hole. 
Here, we explore the effect of the distribution of angular momentum with radius 
in the star, and show that qualitative transitions between high and low angu- 
lar momentum flow, corresponding to high and low luminosity accretion states, 
can effectively be reflected in the energy output, leading to variability and the 
possibility of quiescent times in LGRBs. 



Subject headings: accretion, accretion disks 
namics — supernovae: general 



gamma rays: bursts — hydrody- 



1. Introduction 



The collapse of massive stellar cores producing Supernovae (SN) has been clearly linked 
to the p roduction of Long Gamma Ray Burst s (LGRBs) in the last ten years (see lWoosley fc Bloom 
( 120061 ) ; iGehrels. Ramirez-Ruiz fc Fox! (120091 ) for reviews) , providing tantalizing clues to the 
progenitors and the environments in which they occur jFruchter et al. I2OO6), as well as 



Prochaska et al 



20091 ). A key 



providing fresh glimpses of the high-redshift universe (e.g., 
ingredient in this context is the stellar rotation rate, impacting the energy release and gen- 
eral outcome of the event following the implosion of the Chandrasekhar-mass iron core. 
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Stellar evolution considerations have shown that it is non-trivial to have rapid rotation in 
the inner regions of the star, as mass loss and magnetic fields both conspire to reduce the 
pre-SN angular momentum. Calculations (jSpruitl l2002t iHeger et al.l 120051 ; lYoon &: Langer 
20051 ; IWoosley & Hegerl 12006k hereafter WH06) indicate that the final distributions of spe- 
cific angular momentum may not easily lead to the production of a centrifugally supported 
accretion disk following collapse, with the envelope experiencing essentially radial infall. This 
has generally been believed to be insufficient to power LGRBs since the lack of shocks and 
the associated dissipation are unable to transform gravitational binding energy into radiation 
effectively. A possible evolutionary channel that avoids some of t he pit falls associated with 
the giant phase was proposed in WH06 and by lYoon &: Langerl (120051 ). in which thorough 
mixing on the main sequence in very massive stars reduces mass and angular momentum 
losses. However, a large class of stars remains in which rotation will be substantially slowed 
in late evolutionary stages, thus influencing the outcome of the collapse. 

The low angular momentum regime in collapsars has not been fully explored. In most 
cases, studies of the attending neutrino cooled accretion flows have considered specific ro- 
tation laws that guarantee by a large margin the formation of a centrifugally supported 
disk, either because the angular velocity is assumed to be nearly Keplerian, or because the 
absolute value of the angular momentum given implies a circularization radius much larger 
than the radius of the innermost stable circular orbit (ISCO), which in General Relativity 



2GM RH /c 2 (MacFadven & Wooslev 1999: PoDham. Woosl 


ev & Frver 1999: Heeer et al. 2000: 


Naravan. Piran & Kumarl2001; 


Proea & Beeelmankood: 


Proea. MacFadven. Armitage & Beeelman 


20031: Lee. Ramirez-Ruiz & Pae 


si 20051 


Fuiimoto et al. 20061 Nasrataki et al. 2007). We Dre- 


viouslv studied (Lee & Ramirez- Ruiz 


2006: Lopez-Camara et al. 2009, hereafter Paper I). 



cases at the threshold for disk formation, where a centrifugally supported disk gives way to 
nearly radial inflow due to relativistic effects, finding that even in t he case of slow rotation, 



enough energy may be available throu gh the formation of dwarf disks ( llllarionov fc Beloborodov 



200 ll : IZalamea &: Beloborodov! 120091 ) in near free fall to power LGRBs. However, the distri- 
butions of specific angular momentum we considered were constant in the equatorial plane, 
and smoothly decreasing towards the rotation axis. This is unrealistic, as the specific angu- 
lar momentum increases outwards in the core and envelope, with marked transitions at the 
boundaries between different layers in the star (see, e.g., Figure 2 in WH06). 

Further, as the newborn black hole (BH) accretes infalling matter, both its mass and 
angular momentum increase, raising the threshold value of the critical angular momentum 
required for the formation of a centrifugally supported disk, J crit = 2r g c. The competition 
with the increase of angular momentum in the infalling material determine s the large scale 
properties of the flow. This was recently pointed out and addressed by iJaniuk Proga 
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( 120081 ) in simple form, generically by Kumar. Narayan fc Johnson f l2008al lbh in studying the 
effect of the distr ibution of angular m omentum in the star on the light curve of LGRBs, and 
more recently by lLindner et al.l ( 120101 ) in two-dimensional simulations of collapsar accretion. 



In this study, we explore how different distributions of angular momentum as a function 
of radius, can have important effects on the qualitative properties of the accretion flow, and 
hence on the accretion rate and global energy release, which we quantify through the neutrino 
luminosity, L v . We pay particular attention to the general form and rate of increase of specific 
angular momentum with radius in the star in this respect, and show that state transitions 
may in principle produce observable consequences in LGRBs relevant to variability and 
periods of quiescence. In section § [2] we describe the assumptions and approximations made 
in our calculations and identify the numerical setup. Section § [3] is devoted to presenting 
the flow transitions which result from different angular momentum distributions, and § H] 
to the conditions necessary for them to occur. In § [5] we combine these, given realistic 
angular momentum distributions, to investigate episodic energy release from accretion. Our 
conclusions and prospects for observability are discussed in § [61 



2. Setup and physics 



The evolution of the infalling core and envelope after the formation of the central black 
hole from the stellar Iron core is followed with the same azimuthally symmetric, two di- 
mensional Smooth Particle Hydrodynamics (SPH) code as in Paper I. The computational 
domain covers the region between spherical radii R{ n = 2 x 10 6 cm and R ont = 2 x 10 9 cm, 
without assuming reflection symmetry with respect to the equatorial plane. The equation of 
state is contains contributions to the total pressure P from an ideal gas of a particles and 
free nucleons in nuclear statistical equilibrium (NSE), P ga . s , black-body radiation, P ra .a, rela- 



tivistic e^ 1 pairs of arbitrary degeneracy, P c ± (IBlinnikov. Dunina-Barkovskaya fc Nadyozhin 



19961 ) and neutrinos. Neutronization, and a yariable electron fraction are computed assum- 
ing c harge neutrality and weak equilibrium ( IBeloborodovl 120031 ; iLee. Ramirez-Ruiz fc Page 
20051 ). depending on whether the fluid is optically thin or thick to neutrinos. The neutrino 
emissivities which dominate the cooling arise from e ± capt ure onto free nucleons and e^ 
annih ilation, and are computed fro m the tabu l ated r esults of lLanganke &: Martmez-Pinedo 
(120011 ) and the fitting functions of lltoh et al.l (119961 ) . respectively. A two-stream approx- 
imation including these processes, and c oherent scattering off free nucleons was imple- 



mented to compute the local cooling rat e (jPopham &: Narayanlll995l ; iDi Matteo et al.l 12002 



Janiuk. Yuan. Perna &: Di Matteoll2007l ) . The neutrino luminosity L v is computed as the vol- 
ume integral of the local neutrino emissivity, given the thermodynamical conditions present 
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(see eq.8 in Paper I). The equations of motion also include a full expression for the viscous 
stress tensor ty in azimut hal symmetry. We comput e the coefficient of viscosity with the 



standard a prescription of IShakura fc Sunyaevi (119731 ). namely t/ v = apc s H, where c s is the 



local sound speed, H = c s Q is the pressure scale height, and fl is the local Keplerian orbital 
frequency. We have used a = 0.1 throughout for the simulations described in this paper. 
Given that this form is only applicable if there is a rotating accretion structure, in the usual 
sense of largely centrifugal s upport with a additional press ure corrections, we use a switch 



(similar to that described by iMacFadyen fc Woosley! (Il999l )) so the viscosity only operates 



when rotation dominates over essentially free infall. In practice this is carried out by com- 
paring the local radial and azimuthal components of velocity, and allowing rj v ^ only if 
Vr/vcf, < 1 (the functional form is made continous to avoid spurious transitional behavior 
when rotational support dominates). 

As a generic initial condition we consider the ID pre-supernova models of WH06. 
Specifically we used model 16TI, a rapidly rotating (t> rot = 390km s" 1 ), low metallicity 
(Z = O.O1Z ) WR star of 16M Q with low mass loss (2M total loss up to the pre-supernova 
stage), with an Iron core of Mp e = 1.6M . This particular model did not consider the 
effects of internal magnetic fields. The density, temperature and radial velocity of the stellar 
core was mapped to two dimensions assuming spherical symmetry, while the Iron core was 
condensed onto a point mass at the origin representing the newly for med black hole and 



produ cing a pseudo-Newtonian potential according to the expression of iPaczyhski fc Wiita 



(119801 ) . (j) = — GMbh/ (R — r g ). This reproduces the existence and location of the innermost 
stable circular orbit (ISCO) of the Schwarschild solution in General Relativity, which is the 
leading effect in terms of separating the potential progenitors of GRBs into those which can 
eventually form a centrifugally supported accretion disk and those that cannot. 

The distributions of specific angular momentum, J in pre-supernova stars, as computed 
in evolutionary codes, generally follow a monotonic increase with radius, with sharp tran- 
sitions at the interfaces between shells where the composition changes abruptly. Thus, to 
study the collapse, we separated the angular momentum distribution into radial and polar 
angle components as J = J{R,6) = J(R)Q(9), where R is the spherical radius and 6 is 
the polar angle, and assumed rigid body rotation on shells, with 0(0) = sin 2 6. For the ra- 
dial component, we considered various functional forms: a) linearly increasing with radius, 
J(R) oc R; b) constant, J(R) = Jo with one or more superimposed sharp increases, or spikes 
to mimic the transitions at shell boundaries; c) the distribution given in WH06 for model 
16TI, multiplied by a normalization factor of order unity. 



-5 - 



3. Flow transitions 

3.1. Linear distributions of specific angular momentum. 

In Paper I, we considered simple, constant distributions of J(R) in order to gauge the 
effect of the absolute value of angular momentum and the general properties of the flow 
when compared to previous work, and provide a further guide. When J(R) < J cr it no 
centrifugally supported disk formed, and as the black hole accreted matter and increased its 
mass, so did the value of J cr it, further inhibiting the creation of a disk. The energy released 
came from the equatorial compression of the infalling gas, with peak neutrino luminosities 
L v ~ 10 51 erg s _1 . On the other hand, when J(R) > Jcrit; a hot shocked torus promptly 
formed, but the mass accretion rate onto the black hole, Mbh ~0.1-0.5M© s" 1 , was such 
that the equatorial angular momentum in the flow was always above J crit , and the disk was 
never destroyed. 

Given the general rising trend in the distribution of specific angular momentum in pre- 
SN models mentioned above, we first considered linearly increasing distributions in J(R) 
by writing J(R) = Jo + m(R/R out ) and computing cases with varying J an d m. For 
low (J <C Jcrit) or high (J ^> J C rit) values of the starting point in the distribution, and 
independently of the value of m, the result was akin to that previously obtained, namely, a 
quasi-radial inflow (QRI) or a long-lived accretion disk around the BH, respectively. 

Intermediate cases, however, here when 0.85 < Jo/J cr it — 115, exhibit diverse behavior 
depending on the rate of increase in J(R), given by m. Setting Jo/ Jcrit = 1-05, for m ~ 2, 
a centrifugally supported torus momentarily appears as the gas encounters the centrifugal 
barrier, and is subsequently accreted after a delay of ~. 0.1 s. This is simply a manif estation of 



the well-known "runaway radial" instability ( lAbramowicz. Calvani fc Nobililll983l ). triggered 



by the rapid rise in black hole mass and thus J cr i t . Moreover, the short-lived disk is not able 
to perturb the upstream gas since the inflow is supersonic downstream of the shock front. 
Once the critical rotation rate rises above that of the angular momentum in the envelope, 
the gas is accreted essentially in a free fall time scale and the neutrino luminosity drops by 
about one order of magnitude, to L v ~ 10 51 erg s -1 . There is thus no "memory" in the flow, 
in the sense that once the transitory disk is destroyed, the subsequent material evolves as if 
the disk had never appeared. Figure [1] shows the neutrino luminosity L v for the cases where 
the angular momentum increases either rapidly (blue line), slowly (green line), or moderately 
(red solid line). 



It is thus clear, as explored initially by iJaniuk &: Progal ( 120081 ). that the competition 



between the rise in angular momentum as one moves outward in the star, and the increase 
in BH mass, plays a fundamental role in the establishment of an accretion structure capable 
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Fig. 1. — Neutrino luminosities (scaled to 10 52 erg s -1 ) for simulations with initially linear 
distributions of specific angular momentum, J(R) = Jo + mR/R out . The curves show the 
cases with Jo/J cr it = 1-05, for which the angular momentum increased rapidly (blue, m = 
3, persistent disk); slowly (green, m — 1, no disk); and moderately (red, transitory disk 
correspondent to the critical value m = 2). 

of releasing a significant amount of energy deep in the gravitational potential well. Further- 
more, even a smooth, monotonically increasing distribution of J(R), can produce temporal 
transitions between qualitatively different accretion states, namely, a QRI and a centrifugally 
supported disk. 
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3.2. Rapid transitions in the angular momentum distribution 

The particular feature of strong jumps in J(R) between neighboring shells found in 
stellar evolution calculations led us next to explore how such a qualitative variation in the 
angular momentum distribution can alter the properties of the accretion flow. For this, we 
considered a constant background distribution, with J < J crit , over which two radial shells 
where J(R) reaches 1.5 J CT it are overlaid. Their width is dR ~ 5 x 10 4 cm, and they are 
separated by ~ 2 dR. 

Since the angular momentum of the material initially reaching the BH is below the 
threshold for disk formation, a QRI promptly forms. Once the fluid in the rapidly rotating 
shells reaches the centrifugal barrier, this gives way to a shocked accretion disk, and both the 
accretion rate and neutrino luminosity rapidly rise (see Figure [2]). The delay within which 
the disk forms is essentially the free fall time of the inner shell, t$ oc p ~ x / 2 ~ 0.1 s. In this 
simulated case, the inter-shell region contains enough material to completely overwhelm the 
first disk, leading to a second episode of QRI. This now lasts until the arrival of the second 
shell, with a delay again given by its proper free fall time. Even though the transitions are 
clearly visible in both the accretion rate and the luminosity, some differences clearly stand 
out. First, while Mbh decreases by less than 30%, the luminosity drops by more than one 
order of magnitude. The discrepancy is due to the fundamental morphological change in 
the accretion flow: the QRI has a high mass accretion rate but an extremely low efficiency 
for converting gravitational potential energy into internal energy (it resembles Bondi accre- 
tion in this sense), while the centrifugally supported accretion disk allows for substantial 
energy conversion and dissipation, leading to higher luminosities. The associated accretion 
efficiencies through neutrinos, r) U)!UX = L„/ Mp,uc 2 are in the expected range of (1-10)% for dif- 



feren tially rotating neutrino cooled disks (IDi Matteo et al.ll2002l ; iLee. Ramirez-Ruiz fc Page 



20051 ). and a factor of 10 or more below this for the QRI, where the emission is due to 



equatorial compr ession of the flow to a deg ree which depends on the value of the angular 



momentum used (ILee fc Ramirez- Ruia 120061 . Paper I). Second, there is a slight delay in the 
correlated temporal variations of L v and Mbh- Those in the former come slightly earlier, 
since the spikes in the accretion rate actually correspond to the mass accreted as the slowly 
rotating shells approach the black hole and the differentially rotating disk is destroyed. 

The duration of each interval is dictated in this case by the initial radial position and 
extent of each shell and the corresponding free fall times. Since both shells here have the 
same radial extent and their initial radii are comparable, the associated intervals of activity 
are comparable. Variations in the form and normalization of the rotation rate as a function 
of radius thus impinge upon the duration of active periods and their relative power output. 



We thus conclude that the presence of rapidly rotating shells in a slowly rotating back- 



-8 - 



ground flow can lead to transitions between the QRI, low-luminosity and the disk, high- 
luminosity states, with delays that are correlated with the initial position and width of such 
shells. This is valid as long as the upper and lower limits in the distribution lie below and 
above the critical rotation rate for the existence of a centrifugal barrier. With this in mind, 
we can now address possible effects in GRB central engines in more realistic way, using the 
particular distributions found in WH06. 

4. State transitions and evolutionary pathways 

Before we directly analyze the behavior resulting from the evolutionary model of WH06, 
let us consider under which conditions an infalling shell can destroy a centrifugally supported 
disk. 

Qualitatively, one would expect that the relative mass of successive collapsing regions 
should play a role in determining the outcome of the flow structure. The limits are easy to 
visualize, as a low mass shell with little or no rotation will likely only perturb an existing, 
massive centrifugally supported disk. Conversely, a massive flow in near free fall, due to its 
high ram pressure, will destroy the small pre-existing disk. So where does the transition 
between the two types of flow lie? To address this issue, we performed generic calculations 
(again with model 16TI), and included variations in angular momentum, so that the masses 
of successive shells with high ( Jugh = 3r g c, leading to a first disk) and low (J\ ow = 1.9r g c, 
leading to QRI) angular momentum were and Mqrj, respectively. Their position and 
width were adjusted to explore a range in // = M d /Mq RI from 0.1-10. An outer shell with 
high rotation will lead invariably to the formation of a centrifugally supported disk at late 
times. Note that since we are always assuming rigid body rotation on shells through 0(6), 
it is only the ratio of the masses that is relevant here. 

For the parameters given above, we find that if fi < 1/3, the low-angular momentum 
QRI fully destroys the pre-existing accretion disk on a dynamical time scale. Later, when 
the second shell with J = J high reaches the centrifugal barrier, a new disk is created, which 
persists as long as the inflow has sufficient rotation. This is illustrated in the upper panel of 
Figure |3j where the velocity field shows the flow morphology, and in particular the formation 
of shocks in the disk at various stages of the collapse. For this regime, after an initial QRI 
(t = ti) before the initial shell with J = Jhi g h reaches the inner boundary, a centrifugally 
supported disk is created (t = t 2 ) and subsequently destroyed (t = t 3 ). The associated 
neutrino luminosity shown in the top panel of Figure [3] exhibits a high-state stage between 
two low-state episodes. 
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For jj, > 1/3, the initial disk is able to absorb the impact of the infalling shell and 
survive as a coherent structure. Thus, when the second shell with rapid rotation approaches 
it merely adds to the existing disk activity, and continues its evolution rather unperturbed 
(bottom panel in Figure |3]). The light curve in this case shows continuous emission at the 
high state for t > £2- 

We note here that the absolute values of Jhi g h and J\ ow represent an additional factor 
which can alter the resulting state transitions. Morphologically, the QRI resembles Bondi 
accretion, where technically J\ ow = 0, and it is certainly inefficient in terms of energy release. 
But even a rotation rate below the critical value can have non-negligible consequences when 
combined with a shell slightly above this threshold, and alter the determination of the critical 
mass ratio \i = M d /M QRI required for the transition to occur. In an attempt to quantify this 
effect, we carried out several additional simulations. In the first case, we used \i = 1/3 and 
Jhigh = 3r g c as before, but set J\ ow = 0, thus forcing the more massive shell into strict radial 
infall. The interaction now resulted in the full destruction of the centrifugally supported 
accretion disk. In the second, we changed the mass ratio to \i = 1 and used Jhi g h = 3r g c and 
J\ow = 0. This time, the centrifugally supported disk survived during the entire simulation. 

A simple explanation of this trend can be given as follows. The critical angular momen- 
tum for the appearance of the centrifugal barrier is J cr i t = 2r g c, and we write J\ ow = Ar g c, 
where A is a normalization parameter. Keeping Jhigh = 3r g c as in the above examples, and 
assuming the result after both shells have approached the black hole is a mixed configu- 
ration where the available angular momentum has been distributed into both components, 
the specific angular momentum is now J'. In order for the disk to survive, we must have 
J' > Jcrit, which translates after a few lines of algebra into a condition on the masses as 
Md > (2 — A)Mqri, or > (2 — A). For A = 1.9, as used in our standard case, this implies 
ji > 1/10, roughly consistent with the factor /i = 1/3 found above. If A = 0, as now tested, 
this yields fi > 2 as a condition for disk survival, i.e., it must be more massive in order to 
avoid destruction. In the first test described above, where Ji ow is drastically reduced, the 
outcome is as expected and the disk disappears. In the second, the rise in disk mass is able 
to offset the drop in Ji ow and the centrifugally supported disk persists. The fact that the 
normalization in mass does not strictly follow our simple analytical prescription is clearly due 
to an oversimplification of the dynamics, but it appears nonetheless to capture the essential 
aspect of the problem. 
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5. Critical distributions of angular momentum and episodic energy release 

The distribution of specific angular momentum for representative models taken from 
WH06 is shown in Figure HI They explored a large set of parameters, considering different 
initial zero age main sequence masses, (12 - 35) M , various initial rotation rates, (1 - 14)r g c 
in J(R), initial metallicities, (0.01 - 1)Z Q , mass loss, and the effects of magnetic fields (see 
Tabled]). The resulting distributions span slightly more than two orders of magnitude, and 
all cases clearly show the same qualitative behavior. The general background is that of a 
power law rise, segmented with sharp drops at the interfaces of the various layers in the 
star — for example, the first drop corresponds to the boundary between the iron core and 
the carbon-oxygen (CO) shell. The resulting layers, with varying degrees of rotation and 
mass, must necessarily interact with one another as they collapse onto the central object. 
As before, we chose model 16TI as our fiducial angular momentum distribution, Ji6Ti- 

We are unable to follow the collapse of the outer regions of the core and envelope, 
covering the entire spatial and temporal range required in the present simulations. This 
is mostly due to our choice of location for the inner boundary, with which we are able to 
resolve the transition from a QRI to a rotationally supported disk at a few Schwarzschild 
radii. Nevertheless, the various time scale associated with the arrival of particular shells at 
the centrifugal barrier can be roughly estimated for these models as the free fall time, and the 
previous analysis of high and low-angular momentum shells extended to this regime, which 
we do in what follows. In reality, pressure support within the envelope is non-negligible, so 
the corresponding times and ma ss accretion rates are somewh at longer and lower, respec- 
tively (see also the discussion 



^y 



MacFadven fc Wooslev! (119991 ) related to the latter). The 



calculations by Lindner et al. ( 201ol ) show that the transitions in flow morphology resulti 



mg 



from the structure of the progenitor can indeed be manifested at late times. 

In Figure [5] we show J(R) i6Ti (green dotted line) and J CT it(-R) (black solid line), the 



Table 1: Representative models from WH06. 



Model 


M ZAMS [M ] 


v rot [km s" 1 ] 


z[z Q ] 


M loss [M ] 


B-field 


16TB 


16 


305 


0.01 


0.71 


no 


16TE 


16 


305 


0.01 


4.02 


no 


120N 


12 


400 


0.1 


1.07 


yes 


16TI 


16 


390 


0.01 


2.0 


no 


12TJ 


12 


380 


0.01 


0.46 


yes 


350A 


35 


380 


0.1 


0.6 


yes 


160H 


16 


325 


0.1 


6.82 


yes 
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critical distribution for the formation of a centrifugal barrier assuming all the matter at r < R 
has been accreted by the black hole. It is clear that up to R* ~ 2x 10 9 cm, J(i?)i6Ti < ^crit(-R)- 
Thus for this case, during the first ten seconds (which corresponds to the free fall time 
required for the material departing from R* to reach the black hole), we would have a 
QRI, low-luminosity phase, as the one seen in the upper panel of Figure [3j Given the 
uncertainties in the normalization of J{R) and the spread seen in Figure HJ we now introduce 
a multiplicative scaling factor of order unity / into J(R) i6Ti an d compute the corresponding 
infall time scales. In fact, / may be tuned, for example so that the condition /x = 1/3, 
discussed in § H]is satisfied or not between the initial disk and the first shell with low angular 
momentum (red solid line in Figure with / = 2). Typically, / is of order unity, although 
re-normalizing some models to have \i — 1/3 requires substantial deviations from this value 
(see Table |2}. 

For a normalization leading to a distribution below that with fi = 1/3, a centrifugally 
supported disk created from the material in the shell in region "I" in Figure [5] would be 
destroyed by the QRI produced by the shell labeled "II" , and finally be followed by a second 
stable disk from region "III" . Using the free fall times for the interfaces between regions I, II 
and III we find that after a one-second delay, an initially high-luminosity phase of 1.1 seconds 
would be followed by a quiescent period of ~ 4.2 s before activity resumed. Note that the 
initial black hole with 1.6M is formed from material with R < 2 x 10 8 cm (which includes 
the first large discontinuity in the rotation rate), so no variable energy output emerges from 
this region. 

For vanishing mass of a centrifugally supported disk (p, — >■ 0), the luminosity is only 



that of the QRI at a level L v ~ 10 51 erg s 1 ( lLee fc Ramirez- Ruizll2006l ; lLopez-Camara et al. 



20091 ). As long as0<ii<l/3is satisfied, as in the case shown in Figure^, the low-angular 
momentum shell has greater mass than the centrifugally supported disk and a first high- 
luminosity transient followed by a quiescent period gradually appears, as seen in Figure El 
We note that their respective durations, t& and t q are independent of the assumed strength 
of viscous transport, as they are purely dynamical features. The initial episode lasts until 
the interface of the high and low-angular momentum material, at the edge of the CO shell, 
reaches the centrifugal barrier, t d ~ t s (r co ). The end of the subsequent quiescent interval, t q , 
corresponds to the infall time of the second shell of rapidly rotating gas (above the threshold 
for centrifugal support). The net result is that the sum between td and t q varies slightly as 
a function of \x for a given initial mass configuration (Figure [6]). The neutrino luminosity in 
the high state is L v ~ 10 52 erg s _1 . On the other hand, when the threshold value \i = 1/3 is 
reached, the quiescent interval vanishes as abruptly as the associated flow morphology. For 
larger values, L v can in principle remain uninterruptedly at the high level as long as there 
is mass feeding the disk. While for the particular case shown here the quiescent period lasts 
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for a few seconds, depending on the initial conditions in the progenitor it can in principle be 
substantially longer, matching those observed in GRB events. 



6. Discussion and prospects for observability 



We wish to stress that while we have presented the neutrino luminosity as a measure of 
energy output, it is by no means the only one possible, and should be viewed here as a proxy 
for central engine activity, like the mass accretion rate with which it is very closely correlated. 
One could equally use M (as long as it occurs in state where energy is dissipated efficiently), 
or the power output in magnetic outflows as a measure of the ability to drive relativistic 
outflows. Our numerical scheme is geared towards appropriate handling of thermodynamics 
and the associated neutrino emission, so it is natural to rely on these properties when making 
quantitative statements. 

It is one thing to have a variation in the energy release close to the BH, and another to 
have it manifest itself in the light curve of a LGRB. Now, the minimum time required for 



activity to lead to high energy emission is about 2 s, the crossing time 
originating in the inner regions to pierce the stellar envelope (see, e . g ., 



1999 : lAlov et al. 



200(1 MacFadven et all l200lt lAloy et all 12002 



or a relativistic outflow 



MacFadven fc Woosley 



Ramirez-Ruiz et ajj 



2002 



Zhang et al.ll2003l ). The crossing time is estimated here simply as t cross = R cnv /(c/\/3). In 
Table [2j we present a summary of our results for the models shown in Figure H] and Table HJ 
computed at the threshold where /x = 1/3 and a quiescent period is marginally present. 

After an initial delay (zWay) between the onset of core collapse and the initial rise in 
neutrino luminosity, a first disk is active for an interval td, followed by quiescence and late 
disk creation. Only models with td > t cross , in boldface in Tabled are potentially observable 
as producing a quiescent period in the high-energy light curve under this scheme. Note 
that / is of order unity, but was scaled both to larger and smaller values in order to find a 
condition at the threshold value, fi = 1/3. 

The most obvious limitation of the current study is the range of time scales directly 
modeled. Nevertheless, we believe the current analysis to be of relevance to the generic be- 
havior in collapsing cores. Within this range, the current calculations clearly show that the 
characteristics of the energy output are closely correlated to the distribution of specific angu- 
lar momentum in the progenitor star. The state transitions are abrupt, with the luminosity 
rising or dropping by more than one order of magnitude, and reflect the naturally short (ms) 
time scales in the vicinity of the accretor. The scali ng to larger radii and longer tim e scales 
can lead to quiescent period such as those studied by iRamirez-Ruiz fc Merlonil (120011 ) . which 
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in th is case would be related to d orma nt periods in the central engine ( IRamirez-Ruiz et al. 



20011 ). IRamirez-Ruiz &: Merlonil ( 120011 ) found that a small fraction, ~ 15%, of long GRBs 
exhibit at least one quiescent interval, with about one quarter of these showing two such 
episodes. Furthermore, the durations of the quiescent and subsequent period of activity 
were found to be directly correlated. We believe it is possible to account generically for this 
behavior under the present picture in at least two relevant aspects. 

First, the number of shell-like jumps in the distribution of angular momentum within 
the star is small (one or two), accounting roughly for the number of transitions observed. 
The second point is related to the correlation itself. A short quiescent interval arises when 
the low- angular momentum shell has a narrow radial extent. Its inner limit at T\ is fixed 
by the transition at the edge of the CO core and remains approximately at the same radius, 
independently of the normalization of angular momentum through the parameter /. The 
outer boundary at r2 on the other hand, varies with / (see Figure [5]). Due to the latter and 
to the fact that J cr i t is a monotonically increasing function of radius, higher overall angular 
momentum will lead to: (i) a short quiescent interval, and (ii) a smaller circularization radius 
for the bulk of the innermost matter within the centrifugally supported disk, r circ [J(r2)]. 
Conversely, when the angular momentum normalization is low, the quiescence period will 
be longer, and the circularization radius for the majority of the innermost material will be 
greater. If the evolution of the late-time disk is driven by the viscous transport of angular 
momentum, the corresponding viscous time scale and the duration of the associated accretion 
episode will scale essentially with the initial disk radius. Thus, a longer quiescent period 
will be followed by a longer period of disk activity. The crucial point is that a lower global 
angular momentum normalization will place the bulk of the initial accreting matter capable 
of forming a centrifugally supported disk at a larger circularization radius. There is in 
addition an effective upper limit for this radius if the gas is to release its binding energy 
effectively, because if it is too large, the density and temperature will not be sufficient for 
neutrino cooling to operate. 

As a final point, we may add that just as not all LGRBs exhibit this behavior, likewise it 
is clear from this work t hat not all progenitors are ca pable of producing the state transitions 



presented here. In fact, iPerna fc MacFadyenl (120101 ) have recently argued that the lack of 



flaring activity in most LGRBs on long time scales is a signature of complete mixing in the 
progenitor, which is fully in agreement with the line of argument presented here. Once the 
deposition of energy in the inner regions of the star has produced a relativistic jet capable 
of tra versing the stellar e nvelope, further variability may be reproduced in the overall light 



curve (IZhang et al.ll2003l ). Determining whether this can power precursor activity is another 
matter, requiring the initial episode of accretion to create a low density polar funnel in the 
star, which remains to be addressed. 
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Fig. 3. — Density (in g cm -3 ) and velocity field showing the flow morphology during the 
collapse in the stellar core. The snapshots, at ti < t 2 < t 3 , are taken from simulations of 
model 16TI of WH06 for varying normalization factors in the radial distribution of specific 
angular momentum. Two possible evolutionary pathways are indicated, (a) Top: When 
/i = Md/MQRi < 1/3, the initial QRI (ti) is followed by the formation (t 2 ), destruction (t 3 ) 
and re-emergence of a centrifugal accretion disk, (b) Bottom: When fi > 1/3, once a disk 
of mass Md forms, it persists throughout the evolution as a coherent structure, surviving 
the infall of a low-angular momentum shell of mass Mqrj unperturbed. The transitions are 
reflected in the neutrino luminosity. 
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Fig. 4. — Distributions of specific angular momentum versus radius for representative pre-SN 
evolutionary models, taken from WH06. The initial rise at low radii is followed by a sharp 
drop, secondary increase and fall, and outer increase into the stellar envelope. Note the form 
similarity and the relative normalization differences, which is due to assumed metallicity, 
mass loss and magnetic field effects (see Table [T] for model parameters) . 
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Fig. 5. — Radial distribution of angular momentum for model 16TI of WH06, Jigti{R) 
(green dotted line) and critical curve for the formation of a centrifugal barrier J cr j t (thin 
black solid line) around the central Schwarzschild black hole with M(R). The thick red line 
shows / x J(i?)i6Ti with / = 2. 
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Fig. 6. — Quiescent period, t q (red line) and initial disk activity, td (black line) episodes as 
a function of the high angular momentum disk to low angular momentum shell mass ratio, 
\i = M d /M sh . If fi < 1/3, t q + t d varies slowly, while for // > 1/3 the quiescent episode 
abruptly disappears and t q = 0. The dashed horizontal line approximately indicates the jet 
crossing time, t cr0 ss, through the stellar envelope. In this case t q is few seconds long, but can 
be significantly longer depending on the progenitor's rotation. 
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Table 2: Model evolution time scales. Events with quiescent periods. 



Model 


^dclay 


s 


a t d [s] 


t q [s 


Across .'S 


/ 


16TB 


0.5 


2.3 


9.6 


1.8 


0.6 


16TE 


1.3 


0.9 


3.0 


1.8 


2.0 


120N 


3.5 


0.6 


1.6 


1.5 


1.9 


16TI 


2.3 


1.1 


4.2 


1.0 


2.0 


12TJ 


0.7 


0.7 


4.1 


1.9 


3.9 


350A 


0.5 


2.4 


12.6 


1.5 


80 


160H 


0.4 


2.1 


1.4 


2.0 


39 



"Delay between onset of collapse and initial rise in L v . 



